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We calculate the corrections to the amputated Green's functions of 4-fermion operators, in 
1-loop Lattice Perturbation theory. The novel aspect of our calculations is that they are carried 
out to second order in the lattice spacing, ff(a 2 ). 

We employ the Wilson/clover action for massless fermions (also applicable for the twisted 
mass action in the chiral limit) and the Symanzik improved action for gluons. Our calculations 
have been carried out in a general covariant gauge. Results have been obtained for several popular 
choices of values for the Symanzik coefficients (Plaquette, Tree-level Symanzik, Iwasaki, TILW 
and DBW2 action). 

We pay particular attention to AF = 2 operators, both Parity Conserving and Parity Violating 
(F stands for flavour: S, C, B). We study the mixing pattern of these operators, to 0{a 2 ), using 
the appropriate projectors. Our results for the corresponding renormalization matrices are given 
as a function of a large number of parameters: coupling constant, clover parameter, number of 
colors, lattice spacing, external momentum and gauge parameter. 
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1. Introduction 

A number of flavour-changing processes are currently under study in Lattice simulations. 
Among the most common examples are the decay K — > %% and K°-K° oscillations. From ex- 
perimental evidence, we know that these weak processes violate the CP symmetry. In theory, the 
calculation of the amount of CP violation in K°-K° oscillations requires the knowledge of Bk ■ 

The Kaon Bk parameter is obtained from the AS = 2 weak matrix element: 

(K \O AS = 2 \K°) 

BK ~ ^{K^sy^^sy^)' ' 

where s and d stand for strange and down quarks, and o hS=1 is the effective 4-quark interaction 
renormalized operator, corresponding to the bare operator: 

o A5=2 = (S%d)(S$d) t # = y M (l - 7s)- (1-2) 

The above operator splits into parity-even and parity-odd parts; in standard notation: A5 = 2 = 
®vv+AA ~ ^v1+av Since the above weak process is simulated in the framework of Lattice QCD, 
where Parity is a symmetry, the parity-odd part gives no contribution to the K°-K° matrix element. 
Thus, we conclude that Bk can be extracted from the correlator (xo > 0, yo < 0): 

C K0 K(x,y) = ((^is*)(*)^v^w(0)(£i)ss)(y)>, 0^=^ = (sy M d)(sy M d) + {sy^y 5 d)(sy^y 5 d) , (1.3) 

where * s tne ^ are °P erator an d Oyv+AA * s tne respective renormalized operator. 



In place of the operator in Eq. ( L3 ) it is advantageous to use a four-quark operator with a 



different flavour content (s, d, s' , d'), and with AS = As + As' = 2, namely 0: 

G^T+aa = {sYfid){sy^d') + (sYnY 5 d)(s'rnY5d') + (sy^d'^y^d) + {sy^y 5 d') (sy^y 5 d) , (1.4) 

where now the correlator is given by: Cfc#K>(x,y) = {(dy^s)(x)2Gy^ AA (0)(d'y 5 s'){y)). Making 
use of Wick's theorem one checks the equality: Cx0Kl(x,y) = CKOK(x,y), which means that both 
correlators contain the same physical information. 

The aforementioned matrix elements are very sensitive to various systematic errors. A major 
issue facing Lattice Gauge Theory, since its early days, has been the reduction of effects induced 
by the finiteness of lattice spacing a, in order to better approach the elusive continuum limit. 

In order to obtain reliable non-perturbative estimates of physical quantities (i.e. improving the 
accuracy of Bk) it is essential to keep under control the &(a) systematic errors in simulations or, 
additionally, reduce the lattice artifacts in numerical results. Such a reduction, regarding renormal- 
ization functions, can be achieved by subtracting appropriately the 0{a 2 ) perturbative correction 
terms presented in this paper, from respective non-perturbative results. 

In this paper we calculate the amputated Green's functions and the renormalization matrices 
of the complete basis of 20 four-fermion operators of dimension six which do not need power 
subtractions (i.e. mixing occurs only with other operators of equal dimensions). The calculations 
are carried out up to 1-loop in Lattice Perturbation theory and up to €?(a 2 ) in lattice spacing. Our 
results are immediately applicable to other AF = 2 processes of great phenomenological interest, 
such as D — D or B — B mixing. Let us also mention that in generic new physics models (i.e. beyond 
the standard model), the complete basis of 4-fermion operators contributes to neutral meson mixing 
amplitudes; this is the case for instance of SUSY models (see e.g. [||]). 
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2. Amputated Green's functions of 4-fermion AS = As + As' = 2 operators. 

In this work we evaluate, up to G{a 2 ), the 1-loop matrix element of the 4-fermion operators 1 : 

^ X y = {sXdWYd') = YX E (^M^^X^M^Xw) (2.1) 

0l Y = {sXd')^Yd)=Y^ E (^(x)X klk2 d%( x ))[4 3 (x)Y k3k4 di 4 ^)) (2-2) 

with a generic initial state: d'"*{p^)s'^{p^)\0), and a generic final state: {0\df 2 (p2)sf i l (pi). Spin 
indices are denoted by i, k, and color indices by a, c, d, while X and Y correspond to the following 
set of products of the Dirac matrices: 

X, Y = {1,7 s , Yn,Y^,^ v , 7 s a^^iS^VAJJY <V = 5hi1.IV]- (2.3) 

Our calculations are performed using massless fermions described by the Wilson/clover action. 
By taking nif = 0, our results are identical also for the twisted mass action and the Osterwalder- 
Seiler action in the chiral limit (in the so called twisted mass basis). For gluons we employ a 
3-parameter family of Symanzik improved actions, which comprises all common gluon actions 
(Plaquette, tree-level Symanzik, Iwasaki, DBW2, Liischer-Weisz). Conventions and notations for 
the actions used, as well as algebraic manipulations involving the evaluation of 1-loop Feynman 
diagrams (up to G(a 2 )), are described in detail in Ref. 

To establish notation and normalization, let us first write the tree-level expression for the am- 
putated Green's functions of the operators Gxy and G^ Y : 

A fL(PuP2,P3,P4,r s ,r d ,r s >,r dl )^ 2 2 i ^ 4 a4 =X hh Y hk <5 Ul „ 2 <5 a3 „ 4 , (2.4) 

( AF )fL(PuP2,P3,P4,r s ,r dl r s ,,r d ,)"l" 2 ^ 4 a4 = -X hk Y hh 8 aw 8 a3 a 2 , (2.5) 

where r is the Wilson parameter, one for each flavour. 

We continue with the first quantum corrections. There are twelve 1-loop diagrams that enter 
our 4-fermion calculation, six for each operator Gxy, Gxy- The diagrams d\ —d^ corresponding to 
the operator Gxy are illustrated in Fig. []]. The other six diagrams, d\ — d£, involved in the Green's 
function of G% Y are similar to d\ — d^, and may be obtained from d\ — d(, by interchanging the 
fermionic fields d and d' along with their momenta, color and spin indices, and respective Wilson 
parameters. 

The only diagrams that need to be calculated from first principles are d\ , d2 and d^ , while the 
rest can be expressed in terms of the first three. In particular, the expressions for the amputated 
Green's functions A^J — A^ y can be obtained via the following relations: 

d4 {pUP2,P3,P4,r s ,r d ,r s ,,r d ,)^ i JJ 4 ' i = (A d{ {-p2,-pi,-p4,-p3,r d ,r s ,r d , ,r^) i ^ h \^ i J , (2.6) 
d 5 (PuP2,P3,P4,r s ,r d ,r s r,r d r) h )J 3 f 4 ' i = A dl {p3,p4,P\,P2,r s ',r d ,,r s ,r d )^ i4 \l 2 \ (2.7) 
A d 6 {Pi,P2,P3,P4,r s ,r d ,r sl ,r dl ) h ) 2 j 3 l i ' i = A d} {p3,P4,PhP2,r s > ,r d , ^^r^^J^ - . (2.8) 

Once we have constructed A^J — A^ we can use relation: 

( A )dj (PuP2,P3,P4,r s ,r d ,r s ,,r d r) h )JJ 4 = -A dj <j}i,P4,p3,P2,r s ,r dl ,r sl ,r d ) h \JJ 2 2 , (2.9) 
'The superscript letter F stands for Fierz. 
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S' SJ \ s' s 

P-A /p, p 3 \ /"p, p 3 ' 






sy \ s s 

P, P 3 \ /P, P 3 





P, P 4 }-V^WA P, P« 

t d' A 

di d2 d3 d4 

Figure 1: 1-loop diagrams contributing to the amputated Green's function of the 4-fermi operator XY . 
Wavy (solid) lines represent gluons (fermions). 

to derive the expressions for {A F ) XY — (A F ) X J . From the amputated Green's functions for all twelve 
diagrams we can write down the total 1-loop expressions for the operators &xy and &xy'- 



LA XY 

7=1 



{A F,XY 



1—loop 



E(a f )^ f . 



(2.10) 



In our algebraic expressions for the 1-loop amputated Green's functions A XY , A XY and A XY we 
kept the Wilson parameters for each quark field distinct, that is: r s , rj, r s >, for the quark fields 
s, d, s' and d' respectively. For the required numerical integration of the algebraic expressions of 
the integrands, corresponding to each Feynman diagram, we are forced to choose the square of the 
value for each r parameter. As in all present day simulations, we set: 



2 2 

n = r d 



r 2 - r 2 - 1 



(2-11) 



Concerning the external momenta p\ (shown explicitly in Fig. [T]) we have chosen to evaluate the 
amputated Green's functions at the renormalization point: 



Pl=P2=P3= PA 



(2.12) 



It is easy and not time consuming to repeat the calculations for other choices of Wilson parameters 
and for other renormalization prescriptions. The final 1-loop expressions for A^ y , A XY and A^ y , 
up to &(a 2 ), are obtained as a function of: the coupling constant g, clover parameter c$w-> number 
of colors Nc lattice spacing a, external momentum p and gauge parameter A. 

The crucial point of our calculation is the correct extraction of the full &(a 2 ) dependence 
from loop integrands with strong IR divergences (convergent only beyond 6 dimensions). The 
singularities are isolated using the procedure explained in Ref. [Q]. In order to reduce the number 
of strong IR divergent integrals, appearing in diagram d\ , we have inserted the identity below into 
selected 3-point functions: 



1 / 2 2 _ - 

1 = k + ap +k — ap —2k +162 - sin(£ a ) sin(ap a ) 
ap ^ a 



(2.13) 



where q 2 = 4^^sin 2 (^-) and k(p) is the loop (external) momentum. The common factor in Eq. 
( 2.13 ) can be treated by Taylor expansion. For our calculations it was necessary only to &{a°): 



ap 2 a 2 p 2 



La Pa 

(P 2 ) 2 



(a 2 p 2 ) 



(2.14) 
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Here we present one of the four integrals with strong IR divergences that enter in this calculation: 



K d 4 k sin(k^) sin(£ v 
-n V-X) k 2 k + ap k-ap 



0.002457072288 - '"^^ ^ + a 2 p 2 ^0.00055270353 (6) 



ln(aV) 
~512% 2 ~ 



-a 2 p2 (0.0001282022(1)+^-^) +0.000157122310fl 2 ^^ WpaPv 
M V 768n z J p l i 

\ ( 2 2 \ (p 2 + p 2 ) y 4 

-0.00029731225(4)+ / - 0.000047949674 M „ +0.000268598599 2 " gPf7 



768tt 2 



0.001870841540-^ 



+^(flV) 



the other three integrals can be found in Ref. [|J]. Integrands with simple IR 

•~ handled by well-known techniques. 
3r A^ r and for the special choices: c$w = 0, 



The results for ine oiner mree integrals can De iouna in K.ei. LftJ- iniegranas wiia simpi 
divergences (convergent beyond 4 dimensions) can be handled by well-known techniques. 

Due to lack of space we present only the results for A^ r and for the special choices: c$w 
A = (Landau Gauge), r s = ra = r s > = = 1, and tree-level Symanzik action: 



A^F/ \«|«2"3«4 — S (ft £ _ ° c 'i a 2° 

JX d l \Phih.i 3 k ~ 167 j. 2 I °«1«4°«3"2 



*W'2*W'4 \ fy y 

■2 I "«1«4"«3«2 I X I A 'l'2 1 '3'4 

-0.507914049(6)] + £ (X/V )m , 2 )i 3 ; 4 



I! 



--ln(aV) -0.05294139(2) 

-ln(a 2 /) + 0.0185984988(9) 

O 



0.3977157268533^^ 

P z 



where: 



"(A^i))^ +« 2 (A #(fl 2 ) )J 1 F |, 



(2.15) 



+ l((^7 V )M 2 (^/47 V ) + (^r V )M 2 (W) i 3Hl 
and: 



ip,t (~\a(a 2 p 2 ) +0.09460083(1) 



ip M ^ln(oV) +0.1692905881(6) 



, (2.16) 



(A^(„2) ^ — Xii i 2 Yi 3 i 4 



•} '4 



p 2 f- ^ In (a V) + 1 -32362250(4)^ + 0.06213649(4) ^JL a 
p 2 ^-^ln(aV) +0.059895142(8)^ + 1.01694823(2)/? 2 

PvPji 



E ((Xff)n h Y h u +X h h {Ytf) h , 4 ) x 0.00592406(2) 
li.v 



- E Wn^WU 



Similar expressions exist for A^J and A^ y . 



P M Pv I -^ln(aV) -0.19915360(1) 

p 2 f — ln(aV) - 0.089628048(6) ) - 0.048180850^1 
\240 / pr 

-Pl (- ^ M«V ) + 0. 16608907(3)) 

PuPv I ^ln(aV) -0.21865900(2) +0.140961390— 
\ 360 p 

_ . 1 10138790 ^v + P^v) _ 47763478 ! (8) W 

/? 2 /J 2 



, 2 

-, 



(2.17) 



5 



(a ) corrections to 1-loop matrix elements of 4-fermion operators 



Fotos Stylianou 



3. Mixing and Renormalization of &xy and &Z Y on t ne lattice 



The matrix element (K lOyy+^J-S^) is very sensitive to various systematic errors. The main 
roots of this problem are: a) G(a) systematic errors due to numerical integration, b) the operator 
Oyy+AA mixes with other 4-fermion AS = 2 operators of dimension six. Mixing with operators of 
lower dimensionality is impossible because there is no candidate AS = 2 operator. 

In order to address these problems we have calculated the mixing pattern (renormalization ma- 
trices) of the Parity Conserving and Parity Violating 4-fermion AS = 2 operators (defined below), 
by using the amputated Green's functions obtained in the previous section. A more extensive theo- 
retical background and non-perturbative results, concerning renormalization matrices of 4-fermion 
operators, can be found in Ref. [||] (see also [Q, |6[ [7|]). Next we summarize all important relations 
from Ref. [||] needed for the present calculation. 

One can construct a complete basis of 20 independent operators which have the symmetries 
of the generic QCD Wilson lattice action (Parity P, Charge conjugation C, Flavour exchange sym- 
metry S=(d O d'), Flavour Switching symmetries S'=(s <-^d, / -f-» d') and S"={s <-> d' ,d s')), 
with 4 degenerate quarks. This basis can be decomposed into smaller independent bases according 
to the discrete symmetries P, S, CPS', CPS". Following the notation of Ref. ^ we have 10 Parity 
Conserving operators, Q, (P= + 1 , 5= ±1) and 10 Parity Violating operators, J2, (P= — 1 , S= ±1): 









?VV i 




0Za], \ 




ef 




=\v 


?vv ± 




, 1 


' <bS=±1 


< ef 




=\v 


?ss±l 


? F SS ]-\[0 PP ±6 


ft] , 1 


<aS=±l 






=w 


?ss±t 




ft], 




k 




=\v 




0tt], 







? VA±' 



VA 

ftM[« 



ft] 



; SP'- 



SPi ' 
J SP\ ' 



Summation over all independent Lorentz indices (if any), of the Dirac matrices, is implied. The 
operators shown above are grouped together according to their mixing pattern. This implies that 
the renormalization matrices Z s=±l (3f s=±l ), for the Parity Conserving (Violating) operators, have 
the form: 

/ Z„ Z12 Z13 Z M Z 15 \ S=±1 /iT n \ 5= ' 



7 S=±1 



Z21 Z22 Z23 Z24 Z25 

Z31 Z32 Z33 Z34 Z35 

Z41 Z42 Z43 Z44 Z45 

\Z5i Z52 Z53 Z54 Z55 y 



■S=±l 



/ 2fu 

iT 22 ^23 

^33 

% 4 J°45 

\ J54 iT 55 / 



(3.1) 



Now the renormalized Parity Conserving (Violating) operators, Q s ~ ±l (if S-±1 ), are defined 
via the equations: 

\S=±\ _ ryS=±\ /-vS=±l (?>S=±\ _ av>S=±l G)S=±\ 



0¥>$=±l COA 



(3.2) 



where l,m= 1,...,5 (a sum over m is implied). The renormalized amputated Green's functions 
[s=±\ (j/fS=±i) corresponding to Q s=±l (=S 5=±1 ), are given in terms of their bare counterparts 
L s=±i (jjfS=±i) through: 



fS=±l 



7 -2 7 S=±l T S=±\ 
¥ ^lm ' ' Ll m 1 



S=±l 



?S=±1 



oS=±\ 



(3.3) 
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where Z^ is the quark field renormalization constant. 

In order to calculate the renormalization matrices Z s=±l (3f s=±l ), we make use of the appro- 
priate Parity Conserving (Violating) Projectors p s=±i (£P s=±l ); 



pS=±\ 

p5=±l 
2 

pS=±l 

3 

pS=±l 
4 

d5=±1 



nw+n. 



AA 



64/V c (/V c ±l)' 
U vv 



S=±l 



nvA+n. 



AV 



TIaa _j_ n ss — TIpp 



64(JV?-1) 32/V c (A^-l) 

Ylyy — Haa — IIpp 

32N C (N?-1) 16(/V C 2 -1)' 



32N C (N?-1) ^ 
2W C ±1 

+ ripp 

T 32/V c (/V c 2 -1) 



n 



77 



32N C (N2-1) 

TlfT 



+ 



96N C [N}-1) 
2iV cT l 



a S=±l 



64/V c (/V c ±l)' 

riyA — iIav 



+ 



64(iV2-l) T 


32N C (N?- 


i) 




Hsp — Hps 


32/V c (A£-l) 


16(iV?- 


i) 


n^p + rips 






2iV c ±l 


32/V c (/V c 2 - 


i) 


n^p + rips 







27V cT l 



where n^y = (X^ ® lj 4 i 3 )5 a2ai <5 a4a3 . Again, summation is implied over all independent Lorentz 
indices (if any) of the Dirac matrices. The above Projectors are chosen to obey the following 
orthogonality conditions: 



m (tree) ' 



dl,„, Tr(£?i '-^„(fre e ), 



(3.4) 



where the trace is taken over spin and color indices, and Lt tne y ^ tr< ^ are the tree-level amputated 
Green's functions of the operators Q s=±l , J3 S=±1 respectively. 
We impose the renormalization conditions: 



Sim: Tr(& i ■ ££ m 



S=±l\ 



(3.5) 



By inserting Eqs. (3.3) in the above relations, we obtain the renormalization matrices Z s ±l , 
^es=±\ j n t erms f k nown quantities: 



r S=±l 



where: 



D 



,S=±1 
Im 



Tr(P; 



S=±l 



m / 



(3.6) 



(3.7) 



Note that D s=±1 and ^ s=±1 have the same matrix structure as Z s=±1 and 3f s=±1 respectively. 

Due to lack of space we provide only the matrix @ s=+1 (Parity Violating P = — 1, Flavour 
exchange symmetry 5 = +1) for the special choices: csw = 0, A = (Landau Gauge), r s = = 
r s i =r$ = 1, N c = 3, and tree-level Symanzik action: 



7>S=+1 _ 



/ 1 + 

1+^22 ^23 

^32 1+^33 



V 



1+^44 %5 

^54 1 + ^55/ 



(3.8) 
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where: 

9n = -\ 

S>22 = 1 
^23 = " 
^32 = " 
^33 =4 

^45 = 4 
^4 = 4 
^55 = 4 



167T 2 
167T 2 
167T 2 
167T 2 

167T 2 

jl 
167T 2 

jL 

16k 2 

_l_ 
16k 2 

_l_ 
16k 2 



7.607190(2) +2 ln(flV) + (2.642227(3) - ^ ln(aV)) «V -2.79899088(3) a 2 



{LaP 



2.2995 19(2) + \n(a 2 p 2 ) + f 1 .846794(4) - \- Ha 2 p 2 )] a 2 p 2 - 0.87361421 (2) a - 

V 36 / p l 



2 la Pa 



I. 1931473(4) + (l .4685426(7) ln(aV)) «V - 0.89270364(3) a 
10.970216(2)- 6 ln(a 2 p 2 ) + (6.711307(3) ln(aV)) «V - 1.7027590(1) a 

I I . 595959(2) - 8 ln(aV) + f 3.102499(4) - - ln(a 2 p 2 j) a 2 p 2 + 1 .92846914(2) a 2 ^-^- 

\ 9 / p z 

10.269734(3) -5 ln(aV) - (0.286209(5) - — ln(a 2 p 2 )) a 2 p 2 + 2.01567490(4) a 2 ^^- 

\ 18 / v 



9.732710(2) -5 ln(a 2 /? 2 ) + (4.602710(4) - - ln(fiV)) fl V - 1-07855465(9), 



lLaPo 



1.1783609(6) + - ln(flV) + ( 1.255191(1) - — ln(a 2 p 2 )) a 2 p 2 -0.98220341 (2) a 2 ^^- 
3 V 54 / p L 

iLoPo 



2.297078(2) + y ln(a 2 /? 2 ) + (o.828945(3) - ^ ln(a 2 /)) a 2 /? 2 - 3.06215785(3)fl 



In order to obtain for a given renormalization prescription, one must make use of the inverse 
fermion propagator, S~\ calculated (up to 1-loop and up to G{a ) for massless Wilson/clover 
fermions and Symanzik gluons) in Ref. [Q]. 
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